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The recently measured spin susceptibility of the 
two dimensional electron gas exhibits a strong 
dependence on temperature, which is incompati- 
ble with the standard Fermi liquid phenomenol- 
ogy. Here we show that the observed temper- 
ature behavior is inherent to ballistic two di- 
mensional electrons. Besides the single-particle 
and collective excitations, the thermodynam- 
ics of Fermi liquid systems includes effects of 
the branch-cut singularities originating from the 
edges of the continuum of pairs of quasiparticles. 
As a result of the rescattering induced by in- 
teractions, the branch-cut singularities generate 
non-analyticities in the thermodynamic potential 
which reveal themselves in anomalous tempera- 
ture dependences. Calculation of the spin sus- 
ceptibility in such a situation requires a non- 
perturbative treatment of the interactions. As 
in high-energy physics, a mixture of the collec- 
tive excitations and pairs of quasiparticles can be 
effectively described by a pole in the complex mo- 
mentum plane. This analysis provides a natural 
explanation for the observed temperature depen- 
dence of the spin susceptibility, both in sign and 
magnitude. 

The temperature dependences of the thermodynamic 
quantities in the Fermi liquid have been originally at- 
tributed to the smearing of the quasiparticle distribu- 
tion near the Fermi surface This yields a relatively 
weak, quadratic in temperature, effect. A contribution 
of collective excitations, which in dimensions larger than 
one has a small phase space has been ignored. There 
is a lacuna in this picture. Both the single-particle and 
collective excitations are described by poles in the cor- 
responding correlation functions. However, besides the 
poles there are branch-cut singularities originating from 
the edges of the continuum of pairs of quasiparticles. 
Such branch-cut singularities have not been given ade- 
quate attention in the theory of Fermi liquid systems. In 
the Fermi liquid theory, a rescattering of pairs of quasi- 
particles is considered for the description of the collective 
excitations which exist under certain conditions. This is 
not all that the rescattering of pairs does. Regardless 
of the existence (or absence) of the collective modes, the 
excitations near the edges of the continuum cannot be 
treated as independent as a consequence of the rescat- 
tering. The thermodynamics of Fermi liquid systems is 
not exhausted by the contributions of the single-particle 
and collective excitations. In interacting systems, as a 



result of the multiple rescattering, the branch-cut singu- 
larities generate anomalous temperature dependences in 
the thermodynamic potential. 

Motivated by recent measurements in the silicon 
metal-oxide-semiconductor field-effect transistors (Si- 
MOSFETs) 0, we study here the temperature depen- 
dence of the spin susceptibility, x(T), in the two dimen- 
sional (2D) electron gas in the ballistic regime. Experi- 
ment indicates that in the metallic range of densities and 
for temperatures exceeding the elastic scattering rate, 
T > l/r e j, the electrons in Si-MOSFET behave as an 
isotropic Fermi liquid with moderately strong interac- 
tions. In particular, the Shubnikov-de Haas oscillations 
both without and with an in-plane magnetic field indi- 
cate clearly the existence of a Fermi surface 0, 0, @] . The 
only observation Q incompatible with the simple Fermi 
liquid phenomenology is a surprisingly strong tempera- 
ture dependence of x(^0- This behavior occurs in a wide 
range of densities that rules out proximity to a T = 
quantum critical point as an explanation of the observed 
temperature effect. In this Report we show that such 
a temperature behavior of the spin susceptibility is in- 
herent to 2D ballistic electrons. We explain the experi- 
ment by means of anomalous linear in T terms [|| gener- 
ated by the electron-electron (e-e) interactions in x(T). 
In recent papers, linear in T terms have been studied 
intensely within perturbation theory 0, 0, @, 0, flit - 
However, these works predict the susceptibility increasing 
with temperature, while the trend observed in the exper- 
iment is opposite. Taken seriously, this discrepancy indi- 
cates that we encounter a non-perturbative phenomenon. 
Here we show that a consistent treatment of the effect of 
rescattering of pairs of quasiparticles in different channels 
provides an explanation of the observed temperature de- 
pendence of x- 

How anomalous temperature terms are gener- 
ated in spin susceptibility. Technically, the multiple 
rescattering of pairs of quasiparticles is represented by 
ladder diagrams where each section describes a propa- 
gation of a pair of quasiparticles between the rescatter- 
ing events; see Fig. 1. The collective excitations reveal 
themselves as pole singularities in the ladder diagrams. 
When the pole enters into the continuum of two-particle 
excitations, collective excitation decay. Each of the in- 
termediate sections in the ladder diagrams carries two 
branch-point singularities which reflect the edges of the 
continuum of pairs of quasiparticles. Therefore the cor- 
relation function describing a free propagation of a pair 
of quasiparticles has a branch cut. The analysis of the 



2 



effects of the branch-cut singularities on temperature de- 
pendences in the thermodynamic potential is the object 
of this Report. In the thermodynamic potential the con- 
tribution of the processes of multiple rescattering is given 
by the so-called ring diagrams, i.e., a series of closed lad- 
der diagrams. For the ladder diagrams, the constraints 
imposed by the conservation of the momentum and en- 
ergy are most effective because they are applied to a 
minimal number of quasiparticles. In this way, the dom- 
inant terms are generated in the thermodynamic poten- 
tial. Otherwise summations over a large number of inter- 
mediate states smear out the singularities generated by 
the rescattering processes. 



liquid renormalization of the g-factor. The function So is 
imaginary when (for a given momentum) the frequency 
lies within the continuum of the particle- hole excitations. 
The edges of the particle-hole continuum reveal them- 
selves in So as a branch-point singularities. Since the 
position of the branch cut \lo + A| < v F q depends on 
the magnetic field, the magnetization of the electron gas 
becomes sensitive to the analytical properties of the two- 
particle correlation function near the edge of the con- 
tinuum. The series of ladder diagrams describing the 
rescattering of p-h pairs generates the following term in 
the magnetization (for derivation see Eqs. (12) and (13) 
in the Appendix): 
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FIG. 1: The diagram illustrating two rescattering events. 
Lines connecting the interaction amplitudes F are called a 
"section" . Sections represent a propagation of a pair of quasi- 
particles between rescatterings. Depending on the direction of 
the arrows (omitted here) the section may represent the prop- 
agation of particle-hole (p-h) or particle-particle (Cooper) 
pairs. 



We have to consider series of the ring diagrams in three 
different channels, i.e., in the particle- hole (p-h), the 
particle-particle (Cooper), and the 2/cF-scattering chan- 
nels. The first two channels are standard for Fermi liquid 
theory. The third one is mostly known in connection with 
the Kohn anomaly in the polarization operator [l2T |. We 
start by analyzing the anomalous temperature terms in 
the p-h channel. Within Fermi liquid theory, the e-e in- 
teraction amplitude depends on the angle between the 
incoming and outgoing directions of a scattered parti- 
cle 9 and 8' and is commonly described in terms of the 
angular harmonics. To understand how the anomalous 
temperature terms are generated in the spin susceptibil- 
ity, let us assume for a moment that the zero harmonic, 
To, dominates the interaction amplitude; 7 = Tq. In the 
case of a single harmonic, propagation of a p-h pair is 
described by the angular averaged dynamic correlation 
function*, So = (S(9)) = oj/*J(uj~+ A) 2 - (qv F ) 2 , where 
A = ,g/is(l + Tq)H is the spin split energy induced by an 
external magnetic field H and (1+Fq) describes the Fermi 



We work with the dimensionless static amplitudes known in 
Fermi liquid theory [T3l as T^, while the propagation of a p- 
h pair is described by the dynamic correlation function S(8) = 
lu/(lu + A — qvp cosO); see Eq. 7 in the Appendix. Repulsion 
corresponds to > 0, and Pomcranchuk's instability is at 
r£— 00. 



where Co — lo+A; we temporarily put <7/is(l+ro)/2 equal 
to one. Besides the branch-cut singularities originating 
from the particle-hole continuum, the expression in Eq. 1 
exhibits a pole generated by jlo + \J Co 2 — (qv F ) 2 , which 
determines the spectrum of the collective excitations, i.e., 
the spin- wave excitations [ljj]. Note that the expansion, 
either in 7 or in A, destroys the subtle structure of the de- 
nominator changing its analytical properties. Obviously, 
we encounter a non-perturbative phenomenon. 

In the case of a weak magnetic field, A < T, the 
collective excitations and the continuum of particle-hole 
excitations are not well-separated. Therefore, calcula- 
tions of the thermodynamic quantities, e.g., magneti- 
zation, should be performed with care as the contribu- 
tions from the collective and single-particle excitations 
are not independent. Performing the q- integration by 
contours in the complex g-plane (one should keep in mind 
that the analytical properties in the w-plane differ from 
that in the g-plane), we find that this mixture of ex- 
citations is effectively captured by a pole in the com- 
plex momentum plane. This finding is reminiscent of 
the Regge pole description of the scattering processes in 



high-energy physics 15]. For A / 0, the g-integral is 



non- vanishing only when the pole in the complex g-plane 
(a footprint of the spin-wave excitations) moves into the 
imaginary axis.^ This move occurs within an interval 
—A < lo < — A/(l + 7). At small 7 this interval has 
a width 7A that, by the way, explains why in dSM/dA 
we cannot set A to zero. Only after the (/-integration, we 
get for SM an expression which (at non-zero T) is regular 



t An alternative calculation without referring to the complex q- 
planc is presented in Appendix. 
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both in 7 and A: 



SM = -- 



2e F 



-A/(l+ 7 ) 



(uj + A) coth ^-daj , (2) 



where v is the density of states (per spin) at the Fermi 
surface. Expanding in A, we obtain a linear in T correc- 
tion in the spin susceptibility: 
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T 1 
<$X P -h = -2v— In — — , 

e_F\ 1+7 1 + 7 



(3) 



A comment is in order here. At first glance, a linear in 
T term in x(T) cannot be reconciled with the third law of 
thermodynamics, Sr->o = 0, in view of the Maxwell re- 
lation (dM/dT) H = (dS/dH) T . This observation is the 
core of the statement on the textbook level that the para- 
magnetic behavior with ~ T + const cannot exist at 



sufficiently low temperature; see e.g. Ref. [16J. The well 
known vanishing of the coefficient of thermal expansion 
at the absolute zero has the same origin. In this kind of 
argumentation, it is indirectly assumed that the thermo- 
dynamic potential has a regular expansion in both of its 
arguments around H, T = 0} In fact, Eq. 2 demonstrates 
that the magnetization SM = THm 1 (H /T) has a strong 
dependence on the order of limits H,T — > 0. We see 
from Eq. 2 that SM cx HT when T > H, but for T < H 
the temperature dependence disappears and SM oc H 2 . 
The solution to the conflict with the third law of ther- 
modynamics is that the magnetic field range over which 
SM oc H shrinks to zero as T — > 0, and SM acquires a 
non-linear in H behavior outside this range. At T — > 0, 
which unavoidably brings us into the region T < H, the 
only indisputable condition imposed by the third law is 
limited to vanishing of (dM/dT)n- Evidently, Eq. 2 
complies with this requirement at T/H — > 0. Therefore, 
the existence of a linear in T correction in the spin sus- 
ceptibility is legitimate and may persist down to T — > 0, 
provided that T > H. 

Why spin susceptibility decreases with temper- 
ature. The spin susceptibility as given by Eq. 3 contra- 
dicts the experiment. According to Eq. 3 the spin suscep- 
tibility should increase with T, while in the experiment it 
decreases. Below we offer a resolution to this puzzle. We 
would like first to indicate a subtlety in the thermody- 
namic potential term with two rescattering sections (i.e., 
in the term proportional to T 2 ; see Fig. 2. Obviously in 
the ring diagrams the number of sections is equal to the 
number of the interaction amplitudes). We show below 



t We are not aware of a similar discussion of the thermal expansion 
coefficient (as well as elastic constants) at low temperatures. In 
the context of the spin susceptibility the question has been raised 
by Misawa [l7| who guessed (incorrectly) a non analytic form of 
the thermodynamic potential. 



that the term ~ V 2 in the spin susceptibility is heav- 
ily dominated by the scattering sharply peaked near the 
backward direction, 6 — 9' = it (throughout the paper 
the term "backward scattering" will be used to refer to 
this process). This fact leads to far reaching physical 
consequences, because the diagram with two rescattering 
sections dominated by backward scattering can be read 
in three different ways. Such a diagram can be twisted so 
as to also describe the rescattering in the Cooper chan- 
nel^ or two sections in the 2/cp-scattering channel; see 
Fig. 2 for explanations. Therefore, overlapping of all 
three channels takes place. In order to explain the sign 
of the effect, it is necessary to simultaneously consider 
different channels and to avoid the double counting of 
the contributions generated by different channels. This 
is the central point in our calculation of xCO- 

Before proceeding further, let us outline the conse- 
quences of the overlapping of the three channels which 
takes place on the level of two rescattering sections. In- 
stead of counting the term ~ T 2 in the p-h channel, we 
count it within the Cooper channel where it eventually 
gets killed off by logarithmic renormalizations of the in- 
teraction amplitudes. Therefore, we have to subtract the 
two-section term from S\p-h m Eq. 3 which includes it 
along with higher order terms. As a result of this sub- 
traction Eq. 3 has to be replaced by Eqs. (5) and (6) 
below. In the rest of this section we give the details of 
this procedure. 

Let us first show how the backward scattering arises 
in two p-h rescattering sections. This requires a cal- 
culation of xCO m which the full harmonic content of 
Y(6 — 6') is included. When the amplitudes r„ and 
T m with different harmonics n and m are involved, the 
propagation of the p-h pair between the rescattering 
events is described by the dynamic correlation functions 
Snm = [(w - ^uj 2 -(v F q) 2 )/v F q] ln - ml S Q . Despite the 
nontrivial dependence of S nm on the harmonics indices, 
the contribution to x(T) in the second order of e-e inter- 
action amplitudes acquires a very simple form: 



nm 



n-\-m 



r r 



(4) 



Because J2 n (~ l)Tn is equal to the backward scatter- 
ing amplitude r(7r), this contribution reduces to Sx(2) = 



In fact, an arbitrary number of rescattering sections appears in 
the Cooper channel after such twisting, but only two of them 
are used here for the extraction of the anomalous in temperature 
terms. The role of all other sections is to renormalize logarith- 
mically the e-e interaction amplitudes in the Cooper channel. In 
the text, we refer the term "section" in the Cooper channel only 
to those of them that generate linear in T terms. This allows us 
to speak simultaneously about two sections and the renormalized 
e-e amplitudes without confusion. 
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FIG. 2: The diagrams at the top present the ring diagrams 
in the electron-hole channel. Left diagram shows the two- 
section term which is controlled by the backward scattering; 
the momenta in the four Green's functions are along the same 
direction: 1,2 « +k_F and 3, 4 « — kp. The Green functions 
are numbered to keep track of them after the rearrangement 
in different channels. Right diagram at the bottom shows how 
the two-section term can be read in the 2feF-channel. Here 
the shaded areas represent the interaction amplitudes in the 
2fcF-channel. The lines inside the shaded areas are drawn to 
clarify the spin structure and to indicate the source of the 
relevant renormalizations. Left diagrams at the bottom show 
the result of twisting of the two-section term into the Cooper 
channel. In the series of the Cooper ladder diagrams obtained 
in this way only two sections (marked by numbered Green's 
functions) are responsible for the linear in T term in the spin 
susceptibility. The role of all other sections is to renormalize 
logarithmically the e-e interaction amplitudes. 



(T/£f)vT 2 (it).^ We have checked that exactly the same 
result can be obtained by the calculation of two rescatter- 
ing sections in the Cooper channel, or in the 2/cF-channcl. 
In the calculation of the Cooper channel, we use the an- 
gular harmonics of the particle-particle correlation func- 
tions. Once again, despite the nontrivial dependence of 
these correlation functions on their harmonics indices, 
the result reduces to the backward scattering amplitude; 
details will be published elsewhere [18| . Moreover, this 
calculation yields the same coefficient as in Eq. 4. In the 
case of the 2fcp-channel the presence of r(7r) in Eq. 4 is 
evident, but one has to check the coefficient. On the level 
of two rescattering sections the contributions generated 
in three channels coincide (i.e., in 5\{2) an three channels 
overlap), as we described above. 

We now analyze the problem of the renormalizations of 
the linear in T terms. It is easy to check that unlike the 
case of one-dimensional electrons the higher order 
terms in the 2fcp-scattermg channel are not important 
in 2D. Therefore, we will not discuss this channel fur- 



■ A calculation with the use of angular harmonics has been per- 
formed in [Tl| for the anomalous terms in the specific heat; it 
also leads to the backward scattering amplitude r(-7r). 



ther and concentrate on the interplay between the other 
two channels. Up to this point, the interaction ampli- 
tudes have played a rather passive role in our calcula- 
tions. The peak near the backward scattering direction 
has been generated by the dynamic correlation functions 
describing the propagation of pairs of particles in each 
of the channels. The interaction amplitudes have simply 
supplied a featureless coefficient r(7r) in the two-section 
term. To understand the true role of the e-e interac- 
tion in the anomalous temperature corrections we have 
to abandon the central assumption of the microscopic 
Fermi liquid theory that different sections in the ladder 
diagrams are independent. Indeed, when the rescattering 
is dominated by the backward scattering, a strong de- 
pendence of the interaction amplitude on its arguments 
in the p-h channel emerges from the logarithms in the 
Cooper channel (this is a weak version of the parquet 
known for one-dimensional electrons [l9|). In view of 
this circumstance, in the case of two rescattering sections 
we have to take into consideration the dependence of the 
scattering amplitude T(p, p + q + k, p + k, p + q) on 
the arguments q and k. We resolve the problem of the 
logarithms by moving the term with two rescattering sec- 
tions to the Cooper channel where the logarithmic renor- 
malizations originate. This move is possible because the 
terms with two sections in different channels coincide. 
Note also that as a result of moving the two-section term 
to the Cooper channel we avoid the double counting of 
the two-section term in three different channels. After 
this step, as we discussed earlier, the contribution to the 
spin susceptibility from the p-h channel becomes 

$x' = <5x P -h - 5\{2) ■ (5) 

We now consider the contribution to the spin suscep- 
tibility from the Cooper channel. The rescattering in 
the Cooper channel leads to the logarithmic renormal- 
izations of the interaction amplitudes T%(T) =1^/(1 + 
\nep /T) where are harmonics of the amplitude 
in the Cooper channel. At sufficiently small temper- 
atures, the repulsive amplitudes, > 0, vanish as 
r£(T) w l/ln(e F /T). (We do not consider here the de- 
veloping of the instability for the attractive amplitudes 
(Tc| as it is most likely blocked by the disorder in the sys- 
tem studied in Ref. [2|.) Therefore, the linear in T terms 
generated in the Cooper channel arc suppressed at low 
temperatures. Coming back to the discussion preceding 
Eq. 5, we now see that the logarithmic renormalization 
of the amplitudes in the term with two rescattering sec- 
tions in the p-h channel results in full elimination of this 
term at low enough temperatures. Therefore, for the re- 
pulsive e-e interaction when only zero harmonic is kept, 
the temperature dependence of the spin susceptibility is 
given by 

S X = -2u- ( 7 2 /2 + In -j- + -^-) . (6) 
€f \ 1 + 7 1 + 7/ 
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In contrast to the previous calculations, this result pro- 
vides the sign of the temperature dependence of the spin 
susceptibility which coincides with that observed experi- 
mentally 0]. The expression in Eq. 6 has been obtained 
by summation of the ladder diagrams in two channels and 
taking into consideration the overlap of the two-section 
term. In this way we resolve the puzzle of the sign of the 
temperature trend in x(T)- 

We next note that the intervention of the Cooper 
renormalizations in the p-h channel is effective only for 
the term with two rescattering sections. We have checked 
that the situation with a dominant role of the backward 
scattering is not general and it does not occur for terms 
with more than two rescattering sections. A direct cal- 
culation of the term with three interaction amplitudes 
<5 X( 3) = (T/e F )u J a(9 1 9 2 9 3 )Y(9 1 - 9 2 )T(9 2 - 9 3 )T(9 3 - 
9i)d9id9 2 d9 3l performed with the use of the methods 
sketched above, shows that there is only a weak (loga- 
rithmic) singularity near the backward scattering. This is 
far weaker than the sharp ^-function peak near the back- 
ward direction, 9 — 9' = n, in the case of two rescattering 
sections. It is therefore safe to conclude that, unlike the 
case of the two-section term, the logarithmic renormal- 
izations are ineffective for three and more sections in the 
p-h channel. 

n/x,K 




FIG. 3: The experimental data for n/x expressed in K° and 
depicted as a function of temperature, x ls the spin suscep- 
tibility determined from H = 0.7T spin magnetization; n is 
the density of the 2D electron gas in Si-MOSFET. Numbers 
in the inset are density of the corresponding curves in units of 
10 10 cm -2 . Only a fragment of the data of Ref. Q] correspond- 
ing to the ballistic electrons is presented. [Data courtesy of 
M. Reznikov (Technion, Haifa), used with permission.] 

Relation to experiment. Finally, let us discuss the 
result of our analysis in connection with the measure- 
ment of the spin susceptibility in Si-MOSFET Q. In 
Fig. 3 the data for a quantity n/x are presented, where 
n is the density of the 2D electron gas. We focus here 
on the curves corresponding to the ballistic range of the 



densities. These curves exhibit a noticeable rise with 
temperature at T > 1.5K, which is too strong for the 
conventional Fermi liquid theory; the conventional Fermi 
liquid temperature dependence will be invisible on scales 
used in the plot of Fig. 3. The data correspond to the 
range of densities and temperatures where the transport 
is ballistic. The rising curves in this plot indicate that 
the spin susceptibility decreases with temperature. We 
assume here that this temperature dependence is due to 
the term Sx given by Eq. 6, which we multiply by the fac- 
tor 4 to account for two valleys. At lower temperatures 
the discussed effect of the anomalous temperature cor- 
rections is cut-off by disorder. One can expand n/x with 
respect to the temperature corrections: S(n/x) = T /(t) 
, where f(j) = 7 2 /2-ln(l+7)+7/(l+7). The modifica- 
tion of the spin susceptibility by the Stoner factor (l+To) 
drops out from 5(n/x)- It is cancelled by the two factors 
(l + r ) in Sx ignored so far because in the definition of A 
the combination g/iB(l+ro)/2 has been put equal to one. 
The main advantage of the combination S(n/x) is that its 
temperature dependence is determined by the dimension- 
less interaction amplitudes only. In the discussed range 
of densities, parameter r s is about 3^4 (r s is the ratio of 
the energy of the e-e interaction to the kinetic energy). 
Therefore, one may expect the dimensionless interaction 
amplitudes to have a magnitude ~ 1. Perhaps even a 
few leading harmonics may be involved. For n har- 
monics enter in pairs, T n = r_„, and consequently 7(7) 
should be slightly modified because of mixing between 
r„ and r_„; see Appendix for details. When the ampli- 
tude 7 ~ 1 the function f(j) is of order unity (e.g., for 
7 = 1.5 it is equal to 0.7). The slope of the curves 5 (n/x) 
presented in Fig. 3 is also ~ 1, i.e., of the same order of 
magnitude. Together these facts support our conclusion 
that at low temperature the sign and the magnitude of 
the temperature dependence of the spin susceptibility can 
be explained by the theory of the anomalous corrections 
presented in this Report. At temperatures comparable 
with the Fermi energy the logarithmic suppression of the 
interaction amplitudes in the 5x(2) - term should become 
ineffective. If so, when l/ln(ei?/T) > 7 the temperature 
dependence will change sign leading to a non-monotonic 
spin susceptibility. Unfortunately, the temperature range 
of the existing measurement, T/tp < 0.1, does not allow 
to verify this consequence of our theory. 

To conclude, the thermodynamics of Fermi liquid sys- 
tems is not exhausted by the contributions of the single- 
particle and collective excitations. These two types of 
the excitations are described by the poles in the corre- 
sponding correlation functions. However, the theory of 
Fermi liquid systems is not complete without considera- 
tion of the branch-cut singularities. In interacting sys- 
tems, as a result of the rescattering of quasiparticles, 
the branch-cut singularities generate non-analyticities in 
the thermodynamic potential which reveal themselves in 
anomalous temperature dependences. The observed tern- 
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perature dependence in the spin susceptibility of the 2D 
electron gas can be explained in this way. The mecha- 
nism determining the sign of the anomalous terms in the 
spin susceptibility discussed here may have implications 
for the physics near the quantum critical point at the 
ferromagnetic instability. 



APPENDIX 

Here we present the details of the calculation of the 
anomalous temperature corrections to the spin suscep- 
tibility originating from the particle-hole (p-h) channel. 
The propagation of the pair of the quasiparticlcs with the 
opposite spin projections in the p-h channel is described 
(see § 17 in ref. 13 in the main text) by the two-particle 
correlation function [GG] 9jW ,a = u6e/(u> — 8e), where 
Se = e p + q | — e P x = vj?q— A. Here A = 2h is the relative 
shift of the chemical potential equal to the Zeeman en- 
ergy splitting; h = (gus/tyB. It is convenient to single 
out the dynamic part of this correlation function: 



[GG] q ^,A =-v+ vS(9) q , u ,Al 

S(6) q ,u,,A = — T „ 

uj + A- VFq cos 



(7) 



In this work, the static part of [GG] qiU >,A is absorbed 
in the static Fermi-liquid amplitudes T k and will be not 
considered further. The propagation of a p-h pair is de- 
scribed by the dynamic correlation function S(9). (From 
now on, we will omit vp in the product vpq and consider 
q as measured in energy units.) 

Let us first calculate the contribution from the two 
rescattering sections in the p-h channel to the anomalous 
term in the thermodynamic potential: 

Q 2( A) = -(^)E^/f /^coth^ 

ra,m 

X Im{S n - m (q,LU,A)S 

■m—n (g,w,A)}. (8) 

Because the amplitudes with different harmonics T n and 
T m are involved, the propagation of the p-h pair is de- 
scribed by the angular harmonics of the two-particle cor- 
relation function, S n - m = J ' (d9 /2n)S{0)e l( * n - m ^ , where 



So 



s, 



1^0 



uj — \Ju 2 — q 2 



(9) 

Here we introduce uj — u) + A. Integration over q can be 
done by contours in the complex q plane. This leads to 
the following frequency integral: 



n a (A) = (v/e F )Y, 



Performing the frequency integration for the spin sus- 
ceptibility x = —d 2 il/dh 2 we come to Eq. 4 in the main 
text: 



- / dwur coth— signw. 

16 J 2 6 

(10) 
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We now calculate the ladder diagrams in the p-h chan- 
nel. We start with the zero harmonic; r = 7. For 
completeness we will do it in two different ways. The 
contribution to the thermodynamic potential is equal to 



<Jfi(A) 



(12) 



-(vI2€f) ( duj coth ^— — - Im In — 
w 1 J 2 2tt 1 



1 



+ lSo{q,U, A) 
\vj2tp) I duj coth^— — - Im In ■ x ! 



2 2tt 



juj + ^Juj 2 - q 2 



Here we assume that the frequency uj is slightly shifted 
above the real axes. Eq. 1 in the main text for the mag- 
netization 8M — —dSil/dh follows immediately from this 
expression. To proceed further, we observe that uj and q 
enter in Eq. 1121 only through the combination uj 2 — q 2 . 
Therefore the expression inside the integral vanishes un- 
der the action of uj^d/duj + q^d/dq. We use this to 
write the magnetization M — —dil/dh as an integral of 
the full derivative 



SM = (v/e F ) / — coth — uj 
J 2n 2 



dq 



dq 



-Im In • 



"fUJ + yjuj 2 - q 2 



(13) 



Collecting the contribution at the lower limit of the q 
integral, we obtain 

. r , 1 •> f duj , Buj uj . 

6M = iv Uf) / — coth — cDIm In . 14) 

w 'J 2tt 2 -fuj + uj K ' 



Here we have used that \/ui 2 — q 2 \ q ^o = which corre- 
sponds to the correct analytical structure of the square 
root function in So and Si. For uj slightly above the real 
axis we have 



for 7 > Imln 



(7+1)0; + A 



-7T, -A<C^<- T A- 



0. 



otherwise 



for 7 < Im In ■ 



A 
1+7 



(7 + 1)uj + A [0, otherwise 

(15) 

As a consequence of taking an imaginary part, the fre- 
quency integration is restricted to a narrow frequency 
interval around zero. As a result we obtain Eq. 2 of the 
main text: 

v r A/1+7 0u 

5M = -— (uj + A) coth -—duj . (16) 

2e_F J -a 2 
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It is instructive to reproduce the same result by a more 
powerful (but also more delicate) method of the integra- 
tion by contours in the complex q plane. We return to 
the expression for magnetization 5M: 



SM 



Im 



"qdq 



dui coth 



(3us 



2 

"fUSUS 



2ir (us 2 



'')("fu> + \/Cj 2 -q 2 ) 



(17) 



An important property of this expression is that apart 
from the branch cut on the real axes it also has poles 
in the complex q plane originating from zeros of jus 



\j£o 2 - q" 



Solving the equation jus + \J Co 2 — q 2 = 



we find that the poles never appear inside the branch 
cut in the complex q plane (the branch-cuts are along 
the real axis, covering |q| > w). The poles are either 
somewhere on the real axes between —us and us or appear 
as a pair on the imaginary axis. For 7 > the poles 
exist if —A < us < 0. They are imaginary for — A < 
us < — A/(l + 7). [For 7 < the imaginary poles exist 
for — A/(l + 7) < u) < —A.] We see that the conditions 
that the poles are on the imaginary axes lead to the same 
intervals as in Eq. 1151 above. 

The expression under the integral in Eq. 1171 is an odd 
function of q. This allows us to rewrite the q integral as 
a contour integral in the complex q plane. The contour 
consists of two lines going in the opposite direction above 
and below the real axis. More specifically, the part of the 
contour below the real axis goes in the positive direction 
when us > and in the negative direction when us < 0. 
Thus, 



SM = (i//4ejr) 
qdq 



dus coth 



JUSUS 



Pus 



2m (£2 



q 2 ){^ui + \Jw 2 - q 2 ) 



(18) 



When there are no poles, or they are present but located 
on the real axis q, the integral vanishes. (The contour 
C can be deformed to a big circle where the function 
under the integral drops as ~l/q 2 .) However, the integral 
does not vanish if the poles are on the imaginary axis in 
the complex q plane. This occurs only for the frequency 
intervals discussed above. By deforming the contour and 
taking the residue we reproduce the Eq. 2 in the main 
text. 

So far, only one harmonic in the scattering ampli- 
tude has been iterated within the ring diagrams. More- 
over, for the purpose of clarity it has been assumed 
that zero harmonic amplitude Fo is dominant leading 
to Eq. 6 in the main text. When n =/= 0, harmonics 
enter in pairs because harmonic amplitudes for ±n are 
equal, r„ — r_„. When different harmonics are in- 
volved, the segments representing the iterated harmonics 
7±„ = r n /(l + T n us/ \Jus 2 — q 2 ) have to be connected by 



a correlation function S n — m which represents the "tran- 
sition section" from r„ to T m . For a particular case of 
a pair of harmonics T± n , the transition section is S^n- 
The summation of the ring diagrams for n / can be 
performed by counting how many times the section S^n 
appears 



<5fi(7„,7_„;A) = -(v/2e F ) 

Inl- 



dus coth ^~ 
2 



xlm 



^ln- 

2tt 1 



n 



(19) 



In the following we will limit ourselves to the first mix- 
ing term in Stt in which the transition section S2n appears 
two times; this is reasonable when T n < 1. For clarity, 
we will denote T± n as 7^.. Then, 



-{v/2e F ) 

J+UJ 



duj coth Im 
2 



7_w 



qdq 
2tt 



•in 



j+ls + ^Juj 2 - q 2 7_cj + yfu 



+ v^ 2 



(20) 



The q integral is evaluated by contours in the complex q 
plane 



^A 2 ^-7+7_+^±(A), 
SQ±(A) = -(v/e F ) 1+1 . 



4 

- A/1+7 +, 2 t ,j8a, 

dui us coth — — - — 

—A 2 ( 7 _ 

+ (7+ 7-) 



(21) 
(22) 



7+ 



US + J+UJ 

7+) V^-7+w 



The first term <5f2(A) is a part of the two section term 
which reduces to the square of the backward scattering 
amplitude. It gets killed by the logarithmic renormal- 
izations in the Cooper channel as it has been discussed 
in the main text; we will not keep this term anymore. 
The other term in Eq. 1211 is determined by the poles of 
1 / { r j±ui J r y 'us 2 — q 2 ) when they are on the imaginary axis. 
After passing to a new variable a; = (us/ A + 1)(1+j + )/j + 
we obtain in the limit (3 A <C 1: 

fn±(A) = v tfTw_ ^-^-\ (23) 
2e F (7+ -7-) 

where the function describing the mixing of two harmon- 
ics is 



1+7 Jo 



dx 1 



72; 



x - 1 



1 + 7; + ^ 



(24) 

The last factor reduces the numerical value of 4> n so that 
the correction from the mixing of harmonics is noticeable 
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mostly for first few harmonics with 1. Using the fact 
that the magnitude of amplitudes 7+ and 7_ is equal, we 
obtain 



[3] 



Sx± = -2v(T/e F ) r„— — . 



(25) 



Finally, we obtain for a pair of non-zero harmonics in 



<Wo = -2v(T/e F ) /(r n ) + /(r_„) + r 



#n(r r , 



(26) 

All terms on the right-hand side start with T^. The last 
term describes the modification of the temperature term 
in the spin susceptibility for non-zero harmonics because 
of the mixing of r„ and r_„; compare with Eq. 6 in 
the main text. At T = 1, the value 2f(T = 1) w 0.6, 
whereas for n = 1 ; 2; 10 the corresponding values of 
T 2 (d(f> n {T)/dT) are 0.24; 0.16; 0.04. 

A similar analysis can be performed for the terms with 
three (or more) interaction amplitudes. We check in this 
way that for the higher order terms in the p-h channel the 
scattering does not reduce to the backward scattering. 
Therefore, there is no overlap with the Cooper channel, 
and the intervention of the Cooper channel is irrelevant 
here. 
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